Repurposing FDA-approved drugs to target malaria through inhibition of dihydrofolate reductase in the folate biosynthesis pathway: A prospective approach.

Dihydrofolate reductase (DHFR) is a ubiquitous enzyme that regulates the biosynthesis of tetrahydrofolate among various species of Plasmodium parasite. It is a validated target of the antifolate drug pyrimethamine (Pyr) in Plasmodium falciparum (Pf), but its clinical efficacy has been hampered due to the emergence of drug resistance. This has made the attempt to screen Food & Drug Administration-approved drugs against wild- and mutant PfDHFR by employing an in-silico pipeline to identify potent candidates. The current study has followed a virtual screening approach for identifying potential DHFR inhibitors from DrugBank database, based on a structure similarity search of candidates, followed by absorption, distribution, metabolism, and excretion estimation. The screened drugs were subjected to various parameters like docking, molecular mechanics with generalized born and surface area solvation calculations, and molecular simulations. We have thus identified two potential drug candidates, duloxetine and guanethidine, which can be repurposed to be tested for their efficacy against wild type and drug resistant falciparum malaria.

Oplodiol and nitidine as potential inhibitors of Plasmodium falciparum dihydrofolate reductase: insights from a computational study.

Many natural products have been shown to possess antiplasmodial activities, but their protein targets are unknown. This work employed molecular docking and molecular dynamics simulations to explore the inhibitory activity of some antiplasmodial natural products against wild-type and mutant strains of Plasmodium falciparum dihydrofolate reductase (PfDHFR). From the molecular docking study, 6 ligands preferentially bind at the active site of the DHFR domain with binding energies ranging from -6.4 to -9.5 kcal/mol. Interactions of compounds with MET55 and PHE58 were mostly observed in the molecular docking study. From the molecular dynamics study, the binding of 2 of the ligands-nitidine and oplodiol-was observed to be stable against all tested strains of PfDHFR. The average binding free energy of oplodiol in complex with the various PfDHFR strains was -93.701 kJ/mol whereas that of nitidine was -106.206 kJ/mol. The impressive in silico activities of the 2 compounds suggest they could be considered for development as potential antifolate agents.Communicated by Ramaswamy H. Sarma.

Investigating the anti-cancer potential of pyrimethamine analogues through a modern chemical biology lens.

Pharmacological inhibition of dihydrofolate reductase (DHFR) is an established approach for treating a variety of human diseases, including foreign infections and cancer. However, treatment with classic DHFR inhibitors, such as methotrexate (MTX), are associated with negative side-effects and resistance mechanisms that have prompted the search for alternatives. The DHFR inhibitor pyrimethamine (Pyr) has compelling anti-cancer activity in in vivo models, but lacks potency compared to MTX, thereby requiring higher concentrations to induce therapeutic responses. The purpose of this work was to investigate structural analogues of Pyr to improve its in vitro and cellular activity. A series of 36 Pyr analogues were synthesized and tested in a sequence of in vitro and cell-based assays to monitor their DHFR inhibitory activity, cellular target engagement, and impact on breast cancer cell viability. Ten top compounds were identified, two of which stood out as potential lead candidates, 32 and 34. These functionalized Pyr analogues potently engaged DHFR in cells, at concentrations as low as 1 nM and represent promising DHFR inhibitors that could be further explored as potential anti-cancer agents.

Epistasis and pleiotropy shape biophysical protein subspaces associated with drug resistance.

Protein space is a rich analogy for genotype-phenotype maps, where amino acid sequence is organized into a high-dimensional space that highlights the connectivity between protein variants. It is a useful abstraction for understanding the process of evolution, and for efforts to engineer proteins towards desirable phenotypes. Few mentions of protein space consider how protein phenotypes can be described in terms of their biophysical components, nor do they rigorously interrogate how forces like epistasis-describing the nonlinear interaction between mutations and their phenotypic consequences-manifest across these components. In this study, we deconstruct a low-dimensional protein space of a bacterial enzyme (dihydrofolate reductase; DHFR) into "subspaces" corresponding to a set of kinetic and thermodynamic traits [k_{cat}, K_{M}, K_{i}, and T_{m} (melting temperature)]. We then examine how combinations of three mutations (eight alleles in total) display pleiotropy, or unique effects on individual subspace traits. We examine protein spaces across three orthologous DHFR enzymes (Escherichia coli, Listeria grayi, and Chlamydia muridarum), adding a genotypic context dimension through which epistasis occurs across subspaces. In doing so, we reveal that protein space is a deceptively complex notion, and that future applications to bioengineering should consider how interactions between amino acid substitutions manifest across different phenotypic subspaces.

A rugged yet easily navigable fitness landscape.

Fitness landscape theory predicts that rugged landscapes with multiple peaks impair Darwinian evolution, but experimental evidence is limited. In this study, we used genome editing to map the fitness of >260,000 genotypes of the key metabolic enzyme dihydrofolate reductase in the presence of the antibiotic trimethoprim, which targets this enzyme. The resulting landscape is highly rugged and harbors 514 fitness peaks. However, its highest peaks are accessible to evolving populations via abundant fitness-increasing paths. Different peaks share large basins of attraction that render the outcome of adaptive evolution highly contingent on chance events. Our work shows that ruggedness need not be an obstacle to Darwinian evolution but can reduce its predictability. If true in general, the complexity of optimization problems on realistic landscapes may require reappraisal.

From a binding module to essential catalytic activity: how nature stumbled on a good thing.

Enzymes are complex macromolecules capable of catalyzing a wide variety of chemical reactions with high efficiency. Nonetheless, biological catalysis can be rudimentary. Here, we describe an enzyme that is built from a simple protein fold. This short protein sequence - almost a peptide - belongs to the ancient SH3 family of binding modules. Surprisingly, this binding module catalyzes the specific reduction of dihydrofolate using NADPH as a reducing cofactor, making this a dihydrofolate reductase. Too small to provide all the required binding and catalytic machinery on its own, it homotetramerizes, thus creating a large, central active site environment. Remarkably, none of the active site residues is essential to the catalytic function. Instead, backbone interactions juxtapose the reducing cofactor proximal to the target imine of the folate substrate, and a specific motion of the substrate promotes formation of the transition state. In this feature article, we describe the features that make this small protein a functional enzyme capable of catalyzing a metabolically essential reaction, highlighting the characteristics that make it a model for the evolution of primitive enzymes from binding modules.

Kinetic Barrier to Enzyme Inhibition Is Manipulated by Dynamical Local Interactions in E. coli DHFR.

Dihydrofolate reductase (DHFR) is an important drug target and a highly studied model protein for understanding enzyme dynamics. DHFR's crucial role in folate synthesis renders it an ideal candidate to understand protein function and protein evolution mechanisms. In this study, to understand how a newly proposed DHFR inhibitor, 4'-deoxy methyl trimethoprim (4'-DTMP), alters evolutionary trajectories, we studied interactions that lead to its superior performance over that of trimethoprim (TMP). To elucidate the inhibition mechanism of 4'-DTMP, we first confirmed, both computationally and experimentally, that the relative binding free energy cost for the mutation of TMP and 4'-DTMP is the same, pointing the origin of the characteristic differences to be kinetic rather than thermodynamic. We then employed an interaction-based analysis by focusing first on the active site and then on the whole enzyme. We confirmed that the polar modification in 4'-DTMP induces additional local interactions with the enzyme, particularly, the M20 loop. These changes are propagated to the whole enzyme as shifts in the hydrogen bond networks. To shed light on the allosteric interactions, we support our analysis with network-based community analysis and show that segmentation of the loop domain of inhibitor-bound DHFR must be avoided by a successful inhibitor.

Allosteric regulatory control in dihydrofolate reductase is revealed by dynamic asymmetry.

We investigated the relationship between mutations and dynamics in Escherichia coli dihydrofolate reductase (DHFR) using computational methods. Our study focused on the M20 and FG loops, which are known to be functionally important and affected by mutations distal to the loops. We used molecular dynamics simulations and developed position-specific metrics, including the dynamic flexibility index (DFI) and dynamic coupling index (DCI), to analyze the dynamics of wild-type DHFR and compared our results with existing deep mutational scanning data. Our analysis showed a statistically significant association between DFI and mutational tolerance of the DHFR positions, indicating that DFI can predict functionally beneficial or detrimental substitutions. We also applied an asymmetric version of our DCI metric (DCIasym ) to DHFR and found that certain distal residues control the dynamics of the M20 and FG loops, whereas others are controlled by them. Residues that are suggested to control the M20 and FG loops by our DCIasym metric are evolutionarily nonconserved; mutations at these sites can enhance enzyme activity. On the other hand, residues controlled by the loops are mostly deleterious to function when mutated and are also evolutionary conserved. Our results suggest that dynamics-based metrics can identify residues that explain the relationship between mutation and protein function or can be targeted to rationally engineer enzymes with enhanced activity.

In silico screening of inhibitors against human dihydrofolate reductase to identify potential anticancer compounds.

In all species, dihydrofolate reductase (DHFR) is an essential enzyme that regulates the cellular amount of tetrahydrofolate. Human DHFR (hDHFR) activity inhibition results in tetrahydrofolate depletion and cell death. This property has made hDHFR a therapeutic target for cancer. Methotrexate is a well-known hDHFR inhibitor, but its administration has shown some light to severe adverse effects. Therefore, we aimed to find new potential hDHFR inhibitors using structure-based virtual screening, ADMET prediction, molecular docking, and molecular dynamics simulations. Here, we used the PubChem database to find all compounds with at least 90% structural similarity to known natural DHFR inhibitors. To explore their interaction pattern and estimate their binding affinities, the screened compounds (2023) were subjected to structure-based molecular docking against hDHFR. The fifteen compounds that showed higher binding affinity to the hDHFR than the reference compound (methotrexate) displayed important molecular orientation and interactions with key residues in the enzyme's active site. These compounds were subjected to Lipinski and ADMET prediction. PubChem CIDs: 46886812 and 638190 were identified as putative inhibitors. In addition, molecular dynamics simulations revealed that the binding of compounds (CIDs: 46886812 and 63819) stabilized the hDHFR structure and caused minor conformational changes. Our findings suggest that two compounds (CIDs: 46886812 and 63819) could be promising potential inhibitors of hDHFR in cancer therapy.Communicated by Ramaswamy H. Sarma.

Unveiling the Efficacy of Sesquiterpenes from Marine Sponge Dactylospongia elegans in Inhibiting Dihydrofolate Reductase Using Docking and Molecular Dynamic Studies.

Dihydrofolate reductase (DHFR) is a crucial enzyme that maintains the levels of 5,6,7,8-tetrahydrofolate (THF) required for the biological synthesis of the building blocks of DNA, RNA, and proteins. Over-activation of DHFR results in the progression of multiple pathological conditions such as cancer, bacterial infection, and inflammation. Therefore, DHFR inhibition plays a major role in treating these illnesses. Sesquiterpenes of various types are prime metabolites derived from the marine sponge Dactylospongia elegans and have demonstrated antitumor, anti-inflammation, and antibacterial capacities. Here, we investigated the in silico potential inhibitory effects of 87 D. elegans metabolites on DHFR and predicted their ADMET properties. Compounds were prepared computationally for molecular docking into the selected crystal structure of DHFR (PDB: 1KMV). The docking scores of metabolites 34, 28, and 44 were the highest among this series (gscore values of -12.431, -11.502, and -10.62 kcal/mol, respectively), even above the co-crystallized inhibitor SRI-9662 score (-10.432 kcal/mol). The binding affinity and protein stability of these top three scored compounds were further estimated using molecular dynamic simulation. Compounds 34, 28, and 44 revealed high binding affinity to the enzyme and could be possible leads for DHFR inhibitors; however, further in vitro and in vivo investigations are required to validate their potential.

Exploring the potential and identifying Withania somnifera alkaloids as novel dihydrofolate reductase (DHFR) inhibitors by the AlteQ method.

There is an urgent need to discover and develop novel drugs to combat Mycobacterium tuberculosis, the causative agent of tuberculosis (TB) in humans. Alkaloids have been shown to have wide-ranging therapeutic application and could be ideal candidates for drug development, and research is underway to develop new anti-tubercular drugs from natural sources. In this regard, the current research deals with finding novel lead compounds from the Withania somnifera (WS) plant. Broad health benefits of WS are due to the presence of diverse chemical constituents which include anaferine and anahygrine and which belong to the alkaloid family. In the present study, these two compounds have been theoretically studied to understand their electronic properties using the density functional theory (DFT) at the B3LYP/6-311 + G (d,p) level. HOMO and LUMO properties and molecular electrostatic potential (MEP) surface were calculated. Further, to understand the mechanism of action of these compounds and to identify their putative drug target, molecular docking and dynamics studies were employed against Mycobacterium tuberculosis dihydrofolate reductase (DHFR). It was determined that NADP+ affects stability of the complexes by reducing fluctuations of residues 14-23 and 117-126. It was also found that Ile5 and Gln28 play an important role in complexation. Electron density analysis (using the AlteQ method) of the intermolecular region, analyzing both the anaferin-NADP+ and anahygrin-NADP+ complexes showed that anaferin and anahygrin complexes are more stable in the presence of NADP+. It has been established that in most intermolecular contacts the contribution of the ligand to the electron density is greater than that of NADP+. The present study thus provides an excellent way to analyze the effect of anaferine and anahygrine in essential processes of M. tuberculosis.Communicated by Ramaswamy H. Sarma.

Phosphorylation of Thymidylate Synthase and Dihydrofolate Reductase in Cancer Cells and the Effect of CK2α Silencing.

Our previous research suggests an important regulatory role of CK2-mediated phosphorylation of enzymes involved in the thymidylate biosynthesis cycle, i.e., thymidylate synthase (TS), dihydrofolate reductase (DHFR), and serine hydroxymethyltransferase (SHMT). The aim of this study was to show whether silencing of the CK2α gene affects TS and DHFR expression in A-549 cells. Additionally, we attempted to identify the endogenous kinases that phosphorylate TS and DHFR in CCRF-CEM and A-549 cells. We used immunodetection, immunofluorescence/confocal analyses, reverse transcription-quantitative polymerase chain reaction (RT-qPCR), in-gel kinase assay, and mass spectrometry analysis. Our results demonstrate that silencing of the CK2α gene in lung adenocarcinoma cells significantly increases both TS and DHFR expression and affects their cellular distribution. Additionally, we show for the first time that both TS and DHFR are very likely phosphorylated by endogenous CK2 in two types of cancer cells, i.e., acute lymphoblastic leukaemia and lung adenocarcinoma. Moreover, our studies indicate that DHFR is phosphorylated intracellularly by CK2 to a greater extent in leukaemia cells than in lung adenocarcinoma cells. Interestingly, in-gel kinase assay results indicate that the CK2α' isoform was more active than the CK2α subunit. Our results confirm the previous studies concerning the physiological relevance of CK2-mediated phosphorylation of TS and DHFR.

Arylazopyrazole-Based Photoswitchable Inhibitors Selective for Escherichia coli Dihydrofolate Reductase.

Selective inhibitors of Escherichia coli dihydrofolate reductase (eDHFR) are crucial chemical biology tools that have widespread clinical applications. We developed a set of eDHFR-selective photoswitchable inhibitors by derivatizing the structure of our previously reported methotrexate (MTX) azolog, azoMTX. Substitution of the skeletal p-phenylene group of azoMTX with bulky bis-alkylated arylazopyrazole moieties significantly increased its selectivity toward eDHFR over human DHFR. Owing to the physical properties of arylazopyrazoles, the new ligands exhibited nearly complete Z-to-E photoconversion and high thermostability of Z-isomers. In addition, real-time photoreversible control of eDHFR activity was achieved by alternatively switching the illumination light wavelengths.

A conserved SH3-like fold in diverse putative proteins tetramerizes into an oxidoreductase providing an antimicrobial resistance phenotype.

We present a potential mechanism for emergence of catalytic activity that is essential for survival, from a non-catalytic protein fold. The type B dihydrofolate reductase (DfrB) family of enzymes were first identified in pathogenic bacteria because their dihydrofolate reductase activity is sufficient to provide trimethoprim (TMP) resistance. DfrB enzymes are described as poorly evolved as a result of their unusual structural and kinetic features. No characterized protein shares sequence homology with DfrB enzymes; how they evolved to emerge in the modern resistome is unknown. In this work, we identify DfrB homologues from a database of putative and uncharacterized proteins. These proteins include an SH3-like fold homologous to the DfrB enzymes, embedded in a variety of additional structural domains. By means of functional, structural and biophysical characterization, we demonstrate that these distant homologues and their extracted SH3-like fold can display dihydrofolate reductase activity and confer TMP resistance. We provide evidence of tetrameric assembly and catalytic mechanism analogous to that of DfrB enzymes. These results contribute, to our knowledge, the first insights into a potential evolutionary path taken by this SH3-like fold to emerge in the modern resistome following introduction of TMP. This article is part of the theme issue 'Reactivity and mechanism in chemical and synthetic biology'.

Key interactions of pyrimethamine derivatives specific to wild-type and mutant P. falciparum dihydrofolate reductase based on 3D-QSAR, MD simulations and quantum chemical calculations.

Plasmodium falciparum dihydrofolate reductase-thymidylate synthase (PfDHFR-TS) is an important target enzyme in malarial chemotherapy. An understanding of how novel inhibitors interact with wild-type (wtPfDHFR), quadruple-mutant (qmPfDHFR), and human (hDHFR) enzymes is required for the development of these compounds as antimalarials. This study is focused on a series of des-Cl and m-Cl phenyl analogs of pyrimethamine with various flexible 6-substituents. The interactions of these compounds with DHFR enzymes were investigated by 3 D-QSAR, MD simulations, MM-PBSA, and DFT calculations. CoMFA and CoMSIA models were developed with good predictive abilities for wtPfDHFR and qmPfDHFR. For hDHFR, CoMSIA models combined with clogP descriptor were successfully derived. Binding free energy using MM-PBSA and comparison of per residue decomposition energy analyses with the DFT method at M06-2X/6-31G ++(d,p) level of theory indicated that Asp54 and Phe58 play important roles in the binding of the most potent compound in the series (compound 27) with both wtPfDHFR and qmPfDHFR, whereas Arg59 and Arg122 were additionally found to interact with this inhibitor in qmPfDHFR. For hDHFR, the residues Glu30 and Phe34 but not Arg70, equivalent to Asp54, Phe58, and Arg122 in PfDHFR, also play role in compound 27 binding through strong hydrophobic interactions (Phe34) and hydrogen bond network with Glu30, Ile7, and Val115. From the key interactions identified in the DHFR-inhibitor complexes, a general scheme is proposed for designing new inhibitors selective for PfDHFR that is important for the development of novel antifolate antimalarials.Communicated by Ramaswamy H. Sarma.

Screening of indole derivatives as the potent anticancer agents on dihydrofolate reductase: pharmaco-informatics and molecular dynamics simulation.

Dihydrofolate reductase (DHFR) is a ubiquitous cellular enzyme involved in the biosynthesis of nucleotide and protein precursors, thus, the inhibition of human DHFR can be a promising strategy in cancer treatment. The design of effective anticancer drugs is an urgent need today according to the high spread of cancer. The indole molecule with diverse mechanisms of action and anticancer properties is one of the efficient pharmacophores in drug design. Hence, a virtual library of indole derivatives as a scaffold was selected for designing safer and more effective anticancer drugs against DHFR in this work. All indole derivatives utilized in the library design were selected regarding appreciable tumor growth inhibition. Structure-activity relationship (SAR), docking energy, ADMET (absorption, distribution, metabolism, excretion, and toxicity) parameters, and effective non-covalent interactions were used to identify potential anticancer with indole scaffold. Results showed a higher number of indole moieties provide a strong attachment to the DHFR binding pocket and therefore more effective anticancer activity. The indole scaffold in combination with dichlorobenzene improves DHFR inhibition whereas barbituric acid weakens inhibition activity. In the following to validate the docking results, Molecular dynamics (MD) simulation and molecular mechanics generalized-Born surface area (MM-GBSA) indicated the permanent stability of the selected ligands into the DHFR binding pocket and the key amino acids. Therefore, promising pharmacophores based on indole-DHFR interactions were discovered, and the outcome could be useful in guiding future in vitro and in vivo drug discovery in cancer medicine.Communicated by Ramaswamy H. Sarma.

Phthalide derivatives as dihydrofolate reductase inhibitors for malaria: molecular docking and molecular dynamics studies.

Plasmodium falciparum dihydrofolate reductase enzyme (P. falciparum DHFR) is one of the vital drug targets for malaria treatment, as this protein is indispensable for nucleotide metabolic pathways. This research aimed to discover promising phthalide derivatives against both wild and mutant P. falciparum DHFR enzymes through various computational techniques. The binding affinities were investigated using molecular docking, which showed five compounds having the highest affinity scores against both enzymes compared to the reference compounds. MM-GBSA calculations displayed favourable free binding energy. Moreover, the ADMET properties of the compounds are within acceptable ranges. The stability of the ligand-protein complexes was studied by Molecular Dynamics (MD) simulations. Depending on the results obtained from this research, we propose three compounds to be hit against P. falciparum DHFR activity which could be examined experimentally.Communicated by Ramaswamy Sarma.

Identification and validation of potent inhibitor of Escherichia coli DHFR from MMV pathogen box.

The present study is conducted to find the solution of rising antimicrobial resistance (AMR) in Escherichia coli which is a pathogen responsible for fatal systemic infections in human and animals. The enzyme dihydrofolate reductase (DHFR) is found in all organisms. In this study DHFR of E. coli (ec-DHFR) and human DHFR (h-DHFR) is targeted by novel chemical entities (NCE) from the Pathogen box of Medicines for Malaria Venture, Switzerland (MMV) using molecular modelling. The in-silico studies were further validated by in-vitro assays. The virtual screening of 400 MMV compounds was conducted using PyRx standard tool followed by manual docking of selected compounds by Autodock vina and Ligplot program. The in-silico studies showed good binding energy and strong hydrogen bond in docking of MMV675968 with ec-DHFR and no hydrogen bond with h-DHFR. This was further validated by the Molecular dynamic studies that revealed high binding free energy in ec-DHFR and in-vitro assays that produced good synergy in combination study of MMV675968 with last line (meropenem) and last resort (colistin) antibiotics. The extensive MD simulation and energetic analysis thus concludes that MMV675968 targets ec-DHFR. The combination studies were conducted with MMV675968 and FDA approved drugs against a panel of multidrug resistant Escherichia coli isolates. The synergistic results obtained in combination studies concluded that in-vitro data is consistent with in-silico data and that MMV675968 is a potential lead for future process of antimicrobial drug development against the multidrug resistance E. coli.Communicated by Ramaswamy H. Sarma.

Evaluating the accuracy of the AMBER protein force fields in modeling dihydrofolate reductase structures: misbalance in the conformational arrangements of the flexible loop domains.

Protein flexible loop regions were once thought to be simple linkers between other more functional secondary structural elements. However, as it becomes clearer that these loop domains are critical players in a plethora of biological processes, accurate conformational sampling of 3D loop structures is vital to the advancement of drug design techniques and the overall growth of knowledge surrounding molecular systems. While experimental techniques provide a wealth of structural information, the resolution of flexible loop domains is sometimes low or entirely absent due to their complex and dynamic nature. This highlights an opportunity for de novo structure prediction using in silico methods with molecular dynamics (MDs). This study evaluates some of the AMBER protein force field's (ffs) ability to accurately model dihydrofolate reductase (DHFR) conformations, a protein complex characterized by specific arrangements and interactions of multiple flexible loops whose conformations are determined by the presence or absence of bound ligands and cofactors. Although the AMBER ffs, including ff19SB, studied well model most protein structures with rich secondary structure, results obtained here suggest the inability to significantly sample the expected DHFR loop-loop conformations - of the six distinct protein-ligand systems simulated, a majority lacked consistent stabilization of experimentally derived metrics definitive the three enzyme conformations. Although under-sampling and the chosen ff parameter combinations could be the cause, given past successes with these MD approaches for many protein systems, this suggests a potential misbalance in available ff parameters required to accurately predict the structure of multiple flexible loop regions present in proteins.Communicated by Ramaswamy H. Sarma.

Modeling Ensembles of Enzyme Reaction Pathways with Hi-MSM Reveals the Importance of Accounting for Pathway Diversity.

Conventional quantum mechanical-molecular mechanics (QM/MM) simulation approaches for modeling enzyme reactions often assume that there is one dominant reaction pathway and that this pathway can be sampled starting from an X-ray structure of the enzyme. These assumptions reduce computational cost; however, their validity has not been extensively tested. This is due in part to the lack of a rigorous formalism for integrating disparate pathway information from dynamical QM/MM calculations. Here, we present a way to model ensembles of reaction pathways efficiently using a divide-and-conquer strategy through Hierarchical Markov State Modeling (Hi-MSM). This approach allows information on multiple, distinct pathways to be incorporated into a chemical kinetic model, and it allows us to test these two assumptions. Applying Hi-MSM to the reaction carried out by dihydrofolate reductase (DHFR) we find (i) there are multiple, distinct pathways significantly contributing to the overall flux of the reaction that the conventional approach does not identify and (ii) that the conventional approach does not identify the dominant reaction pathway. Thus, both assumptions underpinning the conventional approach are violated. Since DHFR is a relatively small enzyme, and configuration space scales exponentially with protein size, accounting for multiple reaction pathways is likely to be necessary for most enzymes.